Mesoscopic Fluctuations of Elastic Cotunneling 
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We study mesoscopic fluctuations of the conductance through a quantum dot at the wings of 
the Coulomb blokade peaks. At low temperatures, the main mechanism of conduction is the elastic 
cotunneling. The conductance strongly fluctuates with an applied magnetic field. The magnetic 
correlation field is shown to be controlled by the charging energy, and the correlation function 
has a universal form. The distribution function for the conductance obtained analytically shows a 
non-trivial crossover between the orthogonal and unitary ensembles. 
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Statistical theory of electron eigenfunctions in complex 
sytems was developed in the early fifties Q for the de- 
scription of the absorption spectra of large atoms. Soon 
it was realized that the same random matrix (RMT) ap- 
proach can be applied to the excitation spectrum of small 
metallic grains An important distinction should be 
made between the excitations that preserve the num- 
ber of electrons in the grain, and the excitations that 
change this number. This distinction is most easily seen 
in the case where the introduction of an additional elec- 
tron causes charging of the grain and, thus, is associated 
with the charging energy E c . In this case, the addition 
spectrum shows the gap of the width E c , whereas the 
spectrum of electron-hole excitations is controlled by typ- 
ically much smaller energy scale A (here A is the mean 
level spacing of the grain) . 

The addition spectrum manifests itself in transport 
experiments on electron tunneling between two leads 
through a quantum dot. The electrostatic potential of 
the dot, and therefore the energy of the electron addi- 
tion E, can be controlled by a capacitively coupled gate 
electrode ||. At temperatures T much smaller than E c , 
the presence of the gap results in supression of the con- 
ductance at almost all gate voltages. Only if the gate 
voltage is tuned to one of the discrete points of charge 
degeneracy, i.e. E = 0, is this suppression lifted. This 
phenomenon is known as the Coulomb blockade. Ques- 
tion arises whether and when the results of RMT (which 
does not take into account any charging gap) can be ap- 
plied to the statistics of the conductance in the Coulomb 
blockade regime. Being interesting on its own, this ques- 
tion is also relevant for the recent experimental studies 
of the conductance fluctuations B . The answer to this 
question depends on how close the system is to the charge 
degeneracy point. Right at the degeneracy point, and at 
low temperatures T« A, the transport occurs by means 
of resonant tunneling through a single state. Thus, the 
charging energy E c is irrelevant, and RMT provides an 
adequate description of the distribution function of the 
peak heights |5|(| . Supersymmetry calculation allows 
for a description of the crossover ||, with the magnetic 
field B, between the orthogonal (B — 0) and unitary 



(B — > oo) ensembles. 

Situation changes drastically if the system is tuned 
away from the charge degeneracy point, E > A ~ T. In 
this case the transport is due to the virtual transitions of 
an electron via excited states of the dot (so-called elastic 
cotunneling J^]); many levels with energies exceeding E 
contribute to the tunneling. The superposition of a large 
number of tunneling amplitudes changes the properties 
of the conductance fluctuations, which is studied for the 
first time in this Letter. 

We will show, that the correlation function C(AB) for 
the conductance in the co-tunneling regime is universal 
(i.e. all the dependences for different E can be collapsed 
to a single curve upon rescaling of the magnetic field); 
magnetic correlation field B c is inversely proportional to 
y~E, and thus the charging energy controls the fluctua- 
tions of the elastic cotunneling. Furthermore, the func- 
tional form of C(AB) is entirely different from the known 
results ||. Finally, we will find the distribution function 
of the conductance for all values of the magnetic field. 

The quantum dot attached to two leads is described 
by the Hamiltonian: 



H = H L +H f 



H 7 



(la) 



where the Hamiltonians of the left (L) and right (R) leads 
are given by 



(ib) 



and is the one electron energy measured from the 
Fermi level. Hamiltonian of the dot Hr> has the form 01 



H D = y]$c\c k + E c {h-M) i 



k 



= £4 



{cfc, (1c) 



where £jP describes the one-electron spectrum of the dot, 
and the second term in Ho corresponds to the charging 
energy, and E c = e 2 /2C. Here C is the capacitance of the 
dot, and J\f is the conventional dimcnsionless parameter 
related to the gate voltage V g by N — V g /eC g , with C g 
being the gate capacitance. The tunneling Hamiltonian 
couples the leads with the dot, and it has the form 
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Operators a, b and c in Eqs. (lb)-([ld|) are the correspond- 
ing fermionic operators. 

If tunneling is weak, the charge of the dot n is quan- 
tized. Obviously, degeneracy of the charging energy in 
Eq. ( |Tc| ) corresponds to half-integer values Af m = m + i 
of the dimcnsionlcss gate voltage TV. 

If V g is tuned away from a degeneracy point, it takes a 
finite energy E to add one electron (or hole) to the dot, 



E = E C \N-Af m \, \M — Mm\ < 1/2. 



(2) 



Positive (negative) values of N — N m correspond to the 
electron (hole)-like lowest charged excitations. 

We are considering the strong Coulomb blockade away 
from the resonance. Thus, we employ perturbation the- 
ory in the tunneling Hamiltonian (jld|). The lowest non- 
vanishing contribution to the conductance G is 



G 



2ire 2 
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Amplitude A^ v corresponds to the process in which an 
electron (hole) tunnels from state k in the left lead into a 
virtual state in the dot, and then it tunnels out to state 
p of the second lead. This amplitude is given by 
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The denominator in Eq. (|4|) corresponds to the energy 
of virtual state q involved in the cotunncling process and 
the step-function 6{x) selects the dominating (electron or 
hole) channel. 

In the most realistic case (f§[lO) of point contacts, 
Eqs. (§) and (|) may be further simplified. The tunnel- 
ing matrix elements \t l ^ r q \ 2 do not depend on the indices 
k, q, and can be related to the conductances of the point 
contacts, Gi >r — (2ire 2 /h)vdVi, r \t l ' r \ 2 ] here z^, r ,d are the 
emsemble-averaged densities of states in the leads (I, r) 
and dot (d) respectively. Using these definitions, substi- 

and performing the summa- 
i, we find: 



tuting Eq. (Q) into Eq. 
tion over k and p in Eq. (E 



G — - — ^GiG r |F(R;, R r )| 2 . 



(5) 



The dimensionless function i^R^R,,) contains all the 
information about elastic cotunncling through the dot 
between the point contacts located at R; and R r , 



F(R,R r ) = -]T 



^(R)^(R,0 



E 



[t q {X-Mm)], (6) 



where ip q is the one electron wave function in the closed 
dot. It is useful to rewrite F in terms of the retarded and 
advanced one-electron Green functions G R ' A of the dot, 



F(R,R r ) = — 



duj Ql 



2ni U + E 



[u(M-Af m )}, (7) 



where 



->r,a 
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V^(R r )V; g (R) 
lo — £ q ± iO 



(We put h = 1 in all the intermediate calculations). 

Equations (||) and (^) express the elastic cotunneling 
conductance in terms of the exact electron wavefunctions 
of the dot. These functions vary strongly when the mag- 
netic field is applied to the dot or the shape of the dot 
is changed. Thus, the conductance is a random quantity 
and one should consider different moments of the the 
conductance distribution function. We will employ the 
ensemble averaging, which is equivalent to the averaging 
over applied magnetic field or over the peak index m. 
According to Eqs. (||) and (0), averaged moments of the 
conductance are expressed in the terms of the averaged 
product of the Green functions. It is well known p| , that 
if the dot in the metallic regime (the transport mean free 
path or the size of the dot is much larger than the Fermi 
wavelength), and the relevant energies are much larger 
than A, these products can be related to the generalized 
classical correlators - diffuson V and cooperon C: 



(^, Bl (^)^ B2 (s,r)> = 2n Vd V*l<_% (r,s) 



(8a) 



(SSU ( r > s )Ca (r, «0> = 2W™ (r, s) , (8b) 

where (. . .) stand for the ensemble averaging perfomed 
under the fixed magnetic fields B\,B<2. The averages of 
the type (G R G R ) and (G A G A ) are much smaller and can 
be neglected. If the sample is dirty, so that the motion of 
electrons in the dot is diffusive, the diffuson and cooperon 
(||) satisfy the equations 



V^ B *=5(r-s), (9a) 
C^=8(r-s), (9b) 
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where D is the diffusion constant, and A ± is the vector 
potential due to the magnetic field, V x A* = Bi ± B2. 
For the dot in the ballistic regime, the diffusion operator 
in the l.h.s. of Eqs. (^) should be replaced with the Liou- 
villean operator. Solution of Eqs. (^|) with the condition 
of vanishing normal component of the gauge invariant 
current at the boundary of the dot will enable us to find 
all the relevant correlation functions of the conductance 
and we are turning to this calculation now. 

The averaged cotunneling conductance is obtained im- 
mediately by the averaging of Eq. (||) with the help of 
Eqs. (0) and @. The result is 
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Calculation of the second moment of the conductance is 
performed by taking into account that, because the rel- 
evant energy scale is much larger than A, the average 
of the product of four Green functions can be decoupled 
into a product of the pairwise averaged Green functions, 
which are in turn given by Eqs. (0). This yields 



(G(B)G(B + AB)) 



( GlG r 



(11) 
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The term in the second line of Eq. ( |TT| ) corresponds to the 
square of the average conductance, and the last two terms 
describe the conductance fluctuations. The term in the 
third line of Eq. (|ll]) depends only on AB, [cf. Eq. (9a)] 
and it is present both for the orthogonal (B — 0) and for 
the unitary (B — > oo) ensembles. To the contrary, the 
last term in Eq. (^) dies out for the unitary ensemble. 

It is easily seen from Eq. ( |il| ) that the fluctuations are 
always of the order of the conductance itself. This may 
be understood from the following qualitative considera- 
tion. There are A ~ E/A ^> 1 contributions correspond- 
ing to different eigenstates in the cotunncling amplitude 
([IJ). Assume that the phases of these contributions are 
completely random. Conductance is proportional to the 
modulus squared of the sum of these contributions, and 
thus there are A 2 terms in the conductance. Among 
those, A terms do not fluctuate, and the rest A 2 — A are 
random. These random terms, however, do contribute to 
the fluctuation (SG 2 ), and the number of non- vanishing 
terms in it is A 2 — A. Therefore, the average conduc- 
tance is proportional to A, and its r.m.s. fluctuation is 
~ VA 2 — A ~ A. Thus, conductance in the cotunnel- 
ing regime is not a self-averaging quantity despite a naive 
expectation that a large number of virtual states partici- 
pating in the cotunncling may decrease the fluctuations. 

Equations ( |To| ) and ([ll]) are quite general, i.e. they 
are valid for an arbitrary relation between the energy of 
charged excitation E and Thouless energy Et — fiD / L 2 
(here L is the linear size of the dot). In the most inter- 
esting regime, E < Et, the correlation function of con- 
ductance fluctuations C{AB) acquires a universal form, 
as will be shown below. 

Because E < Et, the characteristic frequency to in 
Eqs. @ is much smaller than the lowest non-zero eigen- 
value of the diffusion operator (which is of the order of 



D/L 2 ). Therefore, only the zero frequency mode can be 
retained in the solutions of Eqs. (|^). This mode cor- 
responds to the probability density homogeneously dis- 
tributed over the dot, and the solution has the form: 
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where S is the area of the dot, <&o = ^irhc/e is the flux 
quantum, and the Thouless energy is given by Et = 
ahD/S, with shape-dependent coefficient a of the order 
of unity. Equation ( |l2|) holds also for ballistic cavities; 
the only difference is that the expression for the Thouless 
energy changes to Et — Ti/rfi, with t// being the time of 
flight of an electron across the dot. Thouless energy can 
be independently measured by studying the correlation 
function of mesoscopic fluctuations for the same dot but 
with the contacts adjusted to the ballistic regime. 

Substitution of Eq. ( p"2| ) into Eq. (|l(]) immediately 
yields the known H result for the averaged conductance 



(G) = 



hGiG r A 
2ne 2 E' 



(13) 



However, the fluctuations 8G(B) = G(B) — (G) are large. 
We find from Eq. (O) with the help of Eqs. (ph: 



(sg{b)6g{b + ab)) jab 

~b7 



A 



2B+AB\ 



B c 



(14) 



where the scaling function A (a;) is given by 



A(x) = 



1 



lnx 2 ln(l + x*) + tt arctana; 2 + iLi 2 (— x*) 



(15) 



with 1^2(21) being the second polylogarithm function 
|]l2| . The asymptoic behavior of function A is A(x) — 
1 + (2x 2 lna; 2 )/7r, for x < 1 and A{x) = {nx 2 )' 2 In 4 x 2 
for x ^> 1. 

The correlation magnetic field B c in Eq. ( |l4| ) is con- 
trolled by the charging energy 



B r = 



$0 

S 



(16) 



It is worth noticing from Eq. (|l3|) and Eq. (|T^) that the 
correlation magnetic field B c drops with approaching a 
charge degeneracy point (in agreement with the recent 
experiment Q), whereas the quantity (G)B 2 remains in- 
variant. This invariance can be easily checked experi- 
mentally. 

Let us present also the expression for the experiman- 
tally measurable correlation function of the conductance 
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fluctuations C(AB) = (SG(B)SG(B + AB))/(6G(B) 2 ). 
For both orthogonal and unitary ensembles we obtain 
from Eq. Q 

C(AB) — A (AB/B C ) , (17) 

function A(x) is defined by Eq. Jl5|). We emphasize that 
the functional form of C(AB) is different from the results 
for the peak heights fluctuations 0, see Fig. |l| 



C(AB) 




FIG. 1. The correlation function C(AB) for the conduc- 
tance fluctuations in the elastic cotunneling regime (solid line) 
and for the peak height fluctuations C = [1 + (AB / B c ) 2 ]~ 2 
(dashed line). For the elastic cotunneling, C(AB) is 
non-analytic at AB — + 0, see the inset, and Eq. (Jjij\). 

As we saw, fluctuations of the conductance are of the 
order of the averaged conductance. Thus the distribu- 
tion function is non-gaussian, i.e. it is not characterised 
by its second moment only. Fortunately, for the elastic 
cotunneling the calculation of all the moments is possible 
which enables us to find the distribution function P(g) 
(here we introduced random variable g = G(B)/(G}). 
Function P(g) is defined as 



<J 
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Calculation of the average (G n (B)) entering into Eq. jlq ) 
is performed in a fashion similar to the derivation of 
Eq. (111]). In the universal regime E < Et, we find 



(G) n 



E 

0<j<n/2 



(n!) 



A^ J 

(j\y(n-2])\ U 



(19) 



where we introduced the short-hand notation A = 
A(B/B C ). Substitution of Eq. (ph into Eq. dl|) yields 



P(g) = 



v 7 ! - A 



exp 
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1 - A 
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(20) 



where Io(x) is the zeroth order modified Bessel function 
of the first kind. In the limiting cases A = 1 and A = 0, 



the distribution function P(g) coincides with the Porter- 
Thomas distribution |l^,|l| for the orthogonal and unitary 
ensembles respectively. 

So far we considered the elastic cotunneling only. It 
dominates over the inelastic processes §] at T < VEA, 
which is the typical regime for the modern experiments 
with semiconductor quantum dots |fl|ic[|. At higher 
temperatures, the main conduction mechanism switches 
to the inelastic cotunneling. Nevertheless, the fluctua- 
tions are still determined by the elastic mechanism for 
(EA) 1/2 < T < (E 2 A) 1/3 . At even higher tempera- 
tures, the inelastic contribution dominates also in the 
fluctuations. Their relative magnitude, however, is small, 
< SG 2 n > I < G in > 2 — A/T. The correlation magnetic 
field is controlled by the temperature rather than by the 
charging energy and therefore is independent on the gate 
voltage. 

In conclusion, we studied the statistics of mesoscopic 
fluctuations of the elastic cotunneling. We showed that 
the correlation magnetic field is controlled by the charg- 
ing energy and the correlation function of the conduc- 
tance is universal. 
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